New metastable Charmonium and the ip' anomaly at CDF 

FE Close 

Rutherford Appleton Laboratory, 
Chilton Didcot, Oxon 0X11 OQX, Great Britain. 
FEC@v2.rl.ac.uk 

RAL-94-093 
August 1994 



Abstract 

Production of metastable charmonium states more massive than the ^(3685) is expected 
in models and could be the source of up to 50% of the t/j observed at large pr at the 
Tevatron. Narrow 2 ^ 2~~cc are predicted at around 3.8GeV and radially excited 2 3 Pi i 2 
may also have suppressed hadronic widths making these states potentially extra sources 
of ?/>(3685). Colour octet components are believed to dominate tp production at the 
Tevatron which suggests that hybrid charmonium production may also be prominent. 
Estimates of hybrid production rates and branching ratios into charmonium suggest that 
metastable hybrids with mass « 4GeV may play an important role in generating the 
observed ^(3685), ^(3095) 
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1 Introduction 



The largest discrepancy between the predictions of the Standard Model and experiment may 
be the recent measurements of charmonium production at large transverse momentum in pp 
collisions at the Tevatron|lJ. The clearest evidence for an anomaly comes from the CDF 
report of a ip' production rate that is a factor of about 30 larger than theoretical expectation, 
too large a discrepancy to be accomodated by tinkering with parameters in the theory and 
suggestive of a source additional to those so far included in the calculations. Production of 
the ip also appears to be enhanced relative to expectations by a factor of 2 or 3 [0; though 
smaller than the ip' enhancement in ratio, this is nonetheless significant and a comparable 
enhancement in absolute magnitude. 

In this letter I bring together recent developments in, at first sight, rather unrelated areas 
of theory and experiment and show how they may provide a consistent picture including a 
possible explanation of the charmonium anomaly. In summary; there are reasons to suspect 
that metastable charmonium states exist more massive than the ip' , that their production 
rates in hadronic interactions are comparable to the x states and that their decays could 
feed the lower mass signals leading to a significant enhancement over the present theoretical 
expectations. Three main candidates occur: 

(i) conventional charmonium states 2 ^,2 , predicted at 3.81 to 3.85GeV g, [| and 
thereby metastable (parity forbids the DD decay), 

(ii) the radial excitation 2 3 Pi(~ ?>.QGeV) whose decay — > DD* is near threshold where 
radial wavefunctions can suppress widths]!], ||, ^| and its partner 2 3 P2(3.9 — A.OGeV) whose 
decays into DD, DD* are similarly suppressed by D-wave phase space and dynamical effects, 

(iii) hybrid charmonium states where the gluonic degrees of freedom are dynamically ex- 
cited in the presence of charmed quarks [f?j|- |p^1 . Some of these states are predicted to 
be metastable if their mass is below DD** threshold (~ 4.3Gey)|], |l3|, [l4j]; recent works 
S 11' 12 1 nave reinforced the predictions @, ||, 12 1 that this mass region is the most likely 



for the manifestation of hybrid charmonium states. 

These three scenarios have competing claims. The first is conservative and almost certain 
to be present at some level: the problem is that the decays into the tp and x charmonia are 
suppressed by phase space and/or wavefunction orthogonality (however, the radiative decay 
of 2 h — > ^(3685) +7 may become detectable if the ^(3685) contains significant ?J D\ or other 
non trivial components in its wavefunction). The radiative transitions 2 3 Pi 5 2 — > 2 3 Si+7 could 
be significant sources of tp(3Q85) if the hadronic width of either 2 3 P2 or 2 3 Pi is suppressed, as 
suggested in some models [||, ||, |l5|.The hybrid excitation on the other hand, is more radical 
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and in consequence more interesting. If such states exist in Nature, then production by hard 
gluons in a hadronic environment favouring heavy flavours, such as at the Tevatron, may be 
rather natural. The pathways to feed ifj and x charmonium states are plentiful and some 
enhancement of ip states may be expected. In any event, the definitive test of all of these 
scenarios is to form invariant mass distributions of ip x an d ip' states in combination with 
7,re7r, 77 in order to isolate any source(s) more massive than t/>(3685). 



Our initial point of departure will be the recent series of papers jig, 17, 18, 1£] on the 
production and decay of heavy quarkonia which have, inter alia, stressed the important role 
that colour-octet components can play, especially in the case of P-wave (and, we conjecture, 
D-wave) quarkonia. The essential physical ideas are as follows. 

Any meson is a superposition of many components 

M = i>QQ + 0{v)^ QQg + ... (1) 

in some of which (e.g.ipQQ g ) the heavy quark and antiquark are in a colour-octet state. The 
probability of the QQg state is of order v 1 (heavy quarks don't easily radiate a gluon;?j 2 ~ ^ 
for charmonium ) and so such components can be neglected for many applications. However, 
as first pointed out in ref [17], these components play an enhanced role in the decay and 



production of P states where the QQg contribution to the quarkonium annihilation competes 
with, and in some cases even dominates, that coming from the QQ. The essential physics is 
that the QQ contribution to annihilations of P-wave quarkonium is suppressed by v 2 owing 
to the angular-momentum barrier which pushes the quarks apart; in the QQg component 
the QQ are in S-wave for which there is no angular momentum suppression though the 
probability to find this component in the wavefunction is only at 0{y 2 ) and hence the QQ 
and QQg contributions are comparable in production or decays. 

These ideas have been extended and applied to charmonium production in P-decays [|l|| 
and to the production of both S'-wave and P-wave charmonia by gluon fragmentation fl§|l . 
As a result of the colour-octet mechanism, the production of x states is predicted to play 
a prominent role, in particular their subsequent decay x ~~ * V* + 7 is predicted to dominate 
over direct production of t/j. There is evidence to support this latter prediction and rather 
general support for the prediction that a major component of heavy quarkonia production at 
high pt > 6GeV is fragmentation (namely the production of high Pt parton, in particular a 



gluon, followed by splitting into a quarkonium state) |20], 21, |2^, |23|| , 

Even though this important colour octet contribution helps to explain in part the large 
J/ifj production recently observed at the Tevatron, it still appears to fall short of explaining 
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the full strength of -0 production, both there ||20| , 21, |22| and also in l?-decay [16|. The 



prediction [16] in the latter case is that BR(B — > + . . .) = 0.23% to be compared with 
(0.74 ± 0.21)% (where the cascade from x states has been accounted for and we have used 
the new PDG values @) 

Production of one or more additional metastable charmonium state (s) with masses greater 
than 3685MeV and with significant branching ratios into the conventional charmonia could 
be responsible for the enhanced signals in the above set of experiments. 

2 Colour Octet States 

Independent of the above, there has been considerable interest in the possible existence of 
hybrid states @, |8], ||] where the gluonic degrees of freedom are excited in the presence of the 
quarks. Early extensive studies || modelled these states as QQg where the quarks were in an 



overall colour octet state; subsequent developments have included lattice simulations [10] and 
modelling in terms of flux tubes connecting the quarks ||. In these latter the conventional 
mesons arise when the flux tube is in its ground state and new hybrid states arise when 
the flux tube is excited. In the case of an unexcited flux tube, the colour can be factored 
onto the QQ such that they are colour singlet; when the colour flux tube is excited such a 
factorisation is non trivial and overlap with the colour octet simulations of refs || may 
arise but a detailed relationship is still unclear. 

In contrast to the conventional quarkonium (eq 1), the hybrid states may be considered 



as g [L9| 



H = iP QQg + 0(vW QQ + ... (2) 

If for the ground state the (QQ)s are in S-wave ||, the production of the hybrid state can 
be dominated by the colour octet mechanisms of refs |l(], [17], ^] in leading order. 

The overall quantum numbers of the hybrids depend on whether the gluon coupling is 
chromo-electric or magnetic ||, forming respectively positive and negative parity families, 
(0, 1, 2) ++ , l" 1 or (0, 1,2) h , 1 . In particular, among the chromomagnetic configurations 
we note the state with J PC = 2 ^ which would be D— wave in colour singlet quarkonium, 
and hence highly suppressed in short distance production (0(v 4 )), but produced in leading 
order in the octet configuration, 

D(g QQ2-+) « i^'(0) 2 « v 4 x R s (0) 2 (3) 
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D(g -> QQgJ~ + ) ~ (2J + l)i?L (4) 



where i?g s is the non-perturbative probability for i n the S-wave colour octet configura- 



tion to fragment into a colour singlet bound state and is expected, following ref [16, 17, 18|, to 
be proportional to the probability of finding the QQg octet component in the hadron wave- 
function. Thus if the QQg "octet modelling" of hybrids|7[ ||, |l9| is a guide, we would expect 
that any 2 ^ charmonium production at the Tevatron will favour the hybrid configuration. 

In the case of chromoelectric gluon couplings, which lead to the same overall quantum 
numbers as those of P-wave charmonium, ref f\M finds for the gluon fragmentation probabil- 
ities 

P^y, fa (2J + 1) / -Rj s (5) 
9 XJ 24m Q 108m Q v ; 

where i?o,i,2 = 5,4, 16 respectively and Hi « 15MeV, H% ks 3MeV. If a s « this implies 
fragmentation probabilities of (0.4, 1.8, 2.4)10 -4 for Xo,l,2 where in each case the octet piece 
dominates. Multiplying the fragmentation probabilities by the appropriate radiative branch- 
ing ratios for \J ~ * ip + 7 OI " (0-7,27, 14)10~ 2 yields the probability of J/ip in a gluon jet of 
approximately 8 x 10 -5 , over an order of magnitude larger than the probability 3 x 10 -6 for 



direct fragmentation of a gluon into J/ip\lS\. 



Insofar as the empirical value H$ m 3MeV subsumes the octet probability 0(v 2 « j), we 
may estimate the probabilities for fragmentation into the hybrid counterparts by rescaling 
the Hi , H$ accordingly 

H hybrid _ \2MeV (6) 
H hybrid _ 3 75Me y ( 7 ) 

in which case the fragmentation probabilities for the hybrids are 



P^ X9 (0 ++ ,l ++ ,2++) = (2,8,13) x lO" 4 (8) 

These give a total probability O(10 -3 ) which would lead to the required enhancement of the 
J/ip signal if B.R.^ip + . . .) « 4 X 10~ 2 . 

A problem is that models put these positive parity hybrid states significantly above the 
charmed hadron pair threshold with consequent strong decay widths O(100)MeV . Therefore 
we find do not anticipate a large (O(10 -2 )) branching ratio into charmonium unless there is 
strong mixing with the nearby radial charmonium states in particular 2 3 Pi t 2(3.9 — A.OGeV) 
(see section 3) 
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By contrast negative parity hybrids could be a copious source of the ip states due to two 
fortunate circumstances common in models. First, models based on constituent gluons^, ||], 



flux tubes|j, 11], quenched heavy quark lattice QCD|To| and QCD sum rules|25|] all anticipate 
negative parity hybrid charmonium states, including (0, 1,2) 1 , typically in the range 
4.2 ± 0.2GeV. In addition the negative parity decays into DD,D*D are suppressed in both 
constituent gluon models p3], Q and in flux tube models § such that the effective strong 
decay threshold is at DD** ~ 4.2 — 4.3GeV (an estimate of their widths leads typically to 
0(lMey)|6|). Thus they offer a priori the possibility of radiative branching ratios into x,*0 
of 0(1CP 3 ) or via hadronic modes at the 10 -2 level if the gluonic tumbledown transitions 



ip(2S) — > 4>(1S) > 50% and predictions of similar rates for ID-states |26|] are any guide. 

Production of these states in 77 processes has been discussed in j0|] but I know of no de- 
tailed calculations analogous to[18| of gluonic fragmentation, factoring out chromomagnetic 



and electric gluon modes. Ref[]l9| found that in 77 decays the colour singlet and octet com- 
ponents compete (for the hybrid decay) and cancel in perturbation theory; in the hadronic 
decays this cancellation does not occur and the colour octet component is dominant. Opti- 
mistically one could take over the results of equation(8) to the negative parity case in which 
case one finds a significant possibility of ip production competing with the existing signal (or 
for the V(3685) a dominant effect). However, we suspect that a more realistic calculation 
will find that the results for the negative parity states will effectively be reduced by 0(v 2 ) 
in probability relative to those for the analogous positive parity states, (the g*g — * QQg 
transition with magnetic J p = 1 + , involves an effective P— wave ( 3 -Pi) coupling relative to 
the 3 Si picture and so we will form our more conservative estimate by reducing the octet 
contribution in equation(5) by 0(v 2 ) « |). 

The singlet contribution in equation(5) also needs to be reconsidered. For the ^ hybrid, 
the colour singlet component is So and so the effective contribution 

IT 1 (0-+)f« JLM^l x0 ( v 2)~27MeV (9) 
2ir nig 

For the case of the exotic 1 ^ and the 2 ^ we need only include the octet piece in our 
estimates (see eq(15) et seq below). So these lead to our "optimistic" and "conservative" 
(0(v 2 ) suppression) estimates for the production probabilities for the negative parity hybrids 
in the QQg picture 

3 x 1(T 4 < P(2~ + ) < 13 x 10" 4 (10) 
2 x 10~ 4 < P(l~ + ) < 4 x 10 -4 (11) 
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1(T 5 < P(0 _+ ) < 1 x 10" 4 (12) 

While these follow rather naturally in the QQg picture they are probably only an order of 
magnitude guide for the flux tube model of hybrids. More detailed study is warranted in the 
latter case as the current estimates are encouraging: in both models the hadronic widths are 
of order lMe^Q due to the effective threshold being that of DD**, and so the feeding of 
-0(3685), ip with significant branching ratios leading to an effective probability 0(10~ 5 — 10~ 4 ) 
may be possible when all states and decay chains are summed. The upper end of this range 
would correspond to this being the dominant contribution to ip production at the Tevatron. 

3 Narrow Charmonium 

In potential models of charmonium, additional narrow states are expected above the 0(3685). 
Most probable are 2 h ,2 predicted at 3.81 to 3.82GeV[Q, ||] and the possibility that 
2 3 Pi i 2(3.9 — A.OGeV) have suppressed hadronic widths due to quantum numbers or nodes 
in form factors manifested in decays near threshold [§J]We consider these in turn. 

An estimate of the production rate for the 2 h follows by assuming this to be dominated 
by the QQg component of the wavefunction. In this case we take the result in eq(10) and 
weigh it by 0(v 2 ) « \, leading to probability 

10~ 4 < P(2~ + ) < 3 x 10" 4 . (13) 

We expect that production of the 2 will be much reduced as there is no obvious octet 
enhancement available, and we shall not consider this state further here. 

The total hadronic width of 2 *~ is expected to be in the range bOkeV — 800/cey. Our 
argument is as follows. 



By analogy with ref[17] we separate the decay into singlet and octet pieces 



,i s a Anal irntctl 
T( 1 D 2 ) = Hf x s - + H 8 x — Ls. (14) 



where |27|| 



H d = ~ 200keV (15) 

2irmq 

leading to a width estimate, in the absence of colour octet contribution, of O(50keV). For 
the octet contribution, we take two values: one is the same as used for the 2 ++ (eq 
5; H 8 = 3.3 ± 0.7MeV) and the other will be to reduce this by 0(v 2 ) w \ for the reasons 
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outlined prior to eq(9). Thus for the octet component alone the width would range between 
120keV < T$ < 750keV. Note that the small (large) widths correlate with small (large) 
production probabilities; hence the ratio of these quantities is less uncertain (and is the 
relevant quantity that will enter in our estimate of the ip yield). 

Combining the singlet and octet contributions leads to a width estimate 

170keV < r(2 -+ ) < 800keV (16) 

(note that this also suggests a hybrid 2 *~ annihilation width of O(lMeV)) 

To estimate the branching ratios into ij)(3685), ip and cascades via 1 D 2 - ► 7 + 1 Pi] -Pi - ► 

ip + . . . we use refs[^7|, ||], rescaled for a 1 D 2 ~ 3.81GeV. 

The 1 D 2 — ► V>(3685) + 7 is suppressed by wavefunction orthogonality unless there is 3 D± 

presence in the ^(3685). The rate is 

T{ 1 D 2 -> V(3685) + 7) = 0{2keV) x Pd (17) 

where pd is the D-wave probability. This implies a branching ratio below 1%. This could 
double the ^(3685) strength but does not seem helpful for explaining a 30- fold enhancement 
unless the intrinsic production rate of the 1 D 2 is larger than our estimates. 

Although this is unlikely to resolve the ^(3685) anomaly it does offer the hope of isolating 
some "missing" charmonium states in the Tevatron and B — ► t/j + . . . data. In this regard, the 



cascades to ^(3095) may be interesting. RefpTfl predicts 



r(>D 2 Jp 1 + 7 ) = m e v( I5 A^ (i 8) 

where k is the photon momentum. The resulting branching ratio could be 50%. According 
to ref||| B.R.( 1 P 1 ^ + . . .) w 1% and so 

B.R.{ 1 D 2 ->V + --0 ~ (4±3)10~ 3 (19) 

These estimates suggest that forming mass plots of ip + ir,2ir,j,r] etc. may isolate x Pi 
and 1 D 2 states, and that these will make some additional contribution to the ip signal at the 
Tevatron but are unlikely alone to explain the full anomaly. 

The final candidate is the production and decay of radial charmonium, in particular 
2 3 Pi j2 -> 2 3 Si + 7. The 2 3 P 2 is predicted to lie 200-300MeV above DD threshold but can 
have a reduced width due to the D-wave phase space and also dynamical effects associated 
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with radial excitations [Q, ||, 15]. The 3 Pi — » also is near threshold and the 5 and D- 

waves are affected by radial wavefunction nodes which can conspire to reduce the width and, 
in consequence, increase the branching ratios 2 3 Pi ; 2 — > 2 3 5i + 7. 

The production probabilities of 2 3 Pj are estimated as for their l 3 Pj counterparts (eq5 
et seq). The magnitude of H\ for the radial is similar to the IP states^]; the magnitude of 
the octet contribution is less certain. If it is reduced by radial effects in proportion to the 
well known colour singlet 5-state wavefunctions there is a partial cancellation with the H\ 
contribution and hence an overall reduction relative to the IP of two( 3 Pi) to three( 3 P2). This 
would imply a probability of 10 -4 for each of these states in gluon jets, which is still a factor of 
100 greater than the direct fragmentation into 0(3685). An order of magnitude increase in the 
■0(3685) signal relative to theory would therefore follow if the PP(2 3 Pi2 — > 2 3 Si +7) ~ 10%. 
(A similar conclusion has been made in refpSfl) 

We can estimate the absolute values of the radiative widths by rescaling the known 



values [24] for the bottomonium system. Allowing a factor of two uncertainty for wavefunction 
differences between cc, bb^ we anticipate T(2 3 P2 C — > 2 3 Si c 7) ~ 80 — 160keV and a similar 
magnitude for the 2 3 Pi case. Thus these states are likely to be serious candidates for the 
■0(3685) enhancement only if their total widths are ~ 1 — 2MeV. 

In flux tube|| or 3 Po models]^] the decays of radially excited states can be artificially 
suppressed if the momenta of the produced DD, DD* coincide with nodes in the form factors. 
The L>-wave decay of 3 P2 and the S, D-wave decays of 3 Pi — ► DD* have widths that are rather 
sensitive to the parameters and masses. If we choose parameters by fitting the hadronic decays 
of ^(3772), 3 L>i and other candidate excited ip states ||l5|] the resulting D-wave contributions 
for the 3 Pj decays are in the range 2 — 5MeV. The .S'-wave contributions are highly sensitive 
to the positions of the wavefunction nodes and hence, in particular, the masses of the states. 
If the mass is taken from refQ, 3 Pi = 3.95GeV then Ts ~ 2.5MeV; if the mass is slightly 
higher (3.96GeV) the S'-wave hits a node and the width is killed. 

Thus model calculations of the hadronic widths suggest that 3 P2 has width ~ 2 — 5MeV, 
and that the 3 Pi width may be similar though this latter result is highly sensitive to the 
choice of parameters. 

Our conclusion is that 1 £ ) 2> 3 D2 charmonium states may be present in the CDF data and 
detectable but are unlikely to explain the 0(3685) enhancement. The radially excited states 
2 3 Pi ; 2 may have suppressed widths enabling them to make a measurable contribution, though 
unlikely on their own to explain the full signal. If the colour octet production mechanism is 



important, as refs|16|, |l7j, |18[ have argued, and if hybrid states containing dominant gluonic 
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excitations exist, then these states are likely to make a significant contribution to the ip 
and ip' signals at the Tevatron, and may also be present in the B — > ip + . . .. While our 
estimates are at this stage only approximate, they do suggest that these various states may 
be present at the right order of magnitude to warrant a serious study. Our results suggest 
that a significant percentage of ip, ip' production at the Tevatron, in l?-decay, or in gluon jets 
at LEP and HERA, could come from higher mass metastable charmonium and that invariant 
mass plots including the ip,ip(3685) may reveal new charmonium states. 
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